Self-Similarity and Energy Dissipation in Stepped Polymer Films 
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The surface of a thin liquid film with nonconstant curvature is unstable, as the Laplace pressure 
drives a flow mediated by viscosity. We present the results of experiments on one of the simplest 
variable curvature surfaces: a stepped polymer film. Height profiles are measured as a function of 
time for a variety of molecular weights. The evolution of the profiles is shown to be self-similar. This 
self-similarity offers a precise measurement of the capillary velocity by comparison with numerical 
solutions of the thin film equation. We also derive a master expression for the time dependence of 
the excess free energy as a function of the material properties and film geometry. The experiment 
and theory are in excellent agreement and indicate the effectiveness of stepped polymer films to 
elucidate nanoscale rheological properties. 



The properties of polymers and indeed all molecules in 
thin films and at interfaces continue to stimulate debate. 
The apparent increased mobility of polymers at interfaces 
in both the liquid fil Q and glassy states [1] , their struc- 
tural properties (3,13 1 as weu as the effects of preparation 
and treatment history are areas of concerted effort. 
While length scales of these systems are frequently just 
a few molecular diameters, the physics governing their 
liquid state evolution is quite general and can be used to 
model, for example, geophysical flows (Toj . 

When a flat liquid film is in contact with a substrate of 
lower surface energy, a hole can nucleate to expose some 
of the underlying surface. This phenomenon induces a 
nonconstant curvature of the film driven by interfacial 
tension and leads to a decrease of the total free e nergy as 
the hole grows. This process is called dewetting [ill [l^. 
There have been many studies concerned with details of 
the shape and evolution of the rim which collects the 
dewetted fluid and much has been learned 

about and from the process of turning a flat film into a 
collection of droplets. 

Several techniques that have been used 20|- 2(1 to pro- 
vide insight into the aforementioned problems !lr9] take 
the opposite approach. They rely on the property that 
curved interfaces have more surface area than flat ones. 
Such curved surfaces can be driven to flatten by the sur- 
face tension, 7, but are mediated by the viscosity, 77. In 
several of the surface tension-driven flow techniques, the 
nanoscale topography is imprinted using a square wave 
pattern that varies in one lateral dimension only ■ 
In these studies, the authors typically consider the evo- 
lution of the amplitude of the perturbation, rather than 
that of the profile shape. As in many of the dewetting 
studies, this Letter is concerned with the details of a 
capillary-driven flow profile, but in a different geometry: 
a stepped film. 



Fig. [lla). They are comparable to the square wave pat- 
terns discussed above [23l426l | in that the height is a func- 
tion of only one spatial dimension and time, h = h(x, t). 
They are different in that they have only one height step. 
This two-dimensional geometry allows us to study the 
broadening in the region over which the height changes 
from hi to h\ + %2 (see Fig. Q]) in isolation, rather than 
having to consider the flow from neighboring steps in 
a periodic geometry. With the use of scaling laws, it 
was shown that by measuring a particular lateral length 
scale, one could measure the relative viscosity [22j . In 
this Letter, we provide a theoretical treatment of the 
governing thin film equation which is in agreement with 
the measured height profiles over a wide range of thick- 
ness combinations and over several orders of magnitude 
in viscosities. We explain in detail the dependencies of 
the energy dissipation on the film geometry and capillary 
velocity, 7/77. 

Films were prepared by spin-coating polystyrene (PS) 
dissolved in toluene onto two types of substrates: 1 x 
1 cm 2 Si wafers (University Wafer) rinsed with ultra- 
pure water (18.2 Mficm, Pall, Cascada LS), methanol, 
and toluene (Fisher Scientific, Optima grade), as well as 
freshly cleaved mica substrates (Ted Pella, Inc.). The 
PS molecular weights were M w = 15.5, 55.5, 118 and 
192 kg/mol with polydispersity indices < 1.07 (Polymer 




In a previous contribution 22} , we explained how to 
prepare the simple geometry of stepped polymer films. 
The samples are prepared as schematically shown in 



FIG. 1: (a) Schematic of the as-prepared samples, (b) After 
annealing above the glass transition temperature, the region 
of transition between the two terrace heights broadens, and a 
flow profile can be measured. 
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FIG. 2: (a) Height as a function of position and time (10 < t < 300 min) for three stepped films with molecular weights as 
indicated. From left to right, {/i.1,/12} = {106, 114}, {89, 89}, and {101,99} nm. The inset shows a detail of the 'dip' region 
for the 118 kg/mol sample, (b) Scaled height as a function of scaled position, according to Eq. ([5}; the inset shows the same 
region as the inset of panel (a) in scaled variables. Whereas the volume constraint uniquely determines the origin, the 192 and 
55.5 kg/mol data have been shifted horizontally and vertically in both (a) and (b) for clarity. 



Source, Inc). The film thicknesses ranged from 30 to 200 
nm. These heights were always much larger than the 
typical size of molecules making up the films. The as- 
cast samples were annealed at 130 °C in a vacuum oven 
(10 -5 mbar) for 24 hours (more than two orders of mag- 
nitude greater than the longest relaxation time of the 
largest polymer used [13] )■ 

The preparation of stepped films proceeded as detailed 
previously |22J. Briefly, PS films with thickness hi on a 
mica substrate were floated onto ultra pure water and 
transferred to a Si wafer. After drying, these wafers 
were split along a crystal axis. The split films were again 
floated onto water and picked up using a Si wafer coated 
with a PS film of height hi, thus creating a bilayer with 
step size hi (see Fig. [IJ . Samples are described as hav- 
ing geometry {hi, hi}. All sections of film observed in 
this study had straight edges as viewed with optical and 
atomic force microscopy (AFM, Veeco Caliber) for dis- 
tances of a minimum of 50 um. Figure Eta) at t = 
illustrates an example of the initial stepped profile; see 
also Ref. 0. 

Prior to measurement of each profile evolution, hi 
was measured (as with all heights reported) using AFM. 
Then, the film was annealed for 10 min at 140 °C on a 
hot stage (Linkam) using a heating rate of 90 °C/min 
in air [28J. After cooling to room temperature (~ 
40 °C/min to below the glass transition, T g ~ 100 °C), 
the step profile was measured. Given h(x, t) as the dis- 
tance between the substrate-polymer and air-polymer in- 
terfaces, h(x,t) — h\ is the quantity measured. For the 



55.5, 118 and 192 kg/mol samples, additional annealing 
for 10 min periods was performed, with the profile evo- 
lution obtained at room temperature. After a final mea- 
surement of the height profile, a scratch in the PS film 
was made down to the substrate to measure hi. 

Figure HJa) shows the time evolution (10 < t < 
300 min) of three stepped polymer films with hi ss hi ~ 
100 nm, each with a different molecular weight. On each 
stepped film, there is a prominent 'dip' on the thin side 
of the film and a 'bump' on the thick side, which is con- 
sistent with our earlier observations [221 ] and can also 
be seen in the long wave limit of the data presented by 
Rognin et al. 25j. The evolution proceeds such that po- 
sitions of the profile extrema are separated over wider 
distances: the film flattens with time. 

The flow profile can be understood from the Laplace 
pressure, which arises due to curvature at the fluid in- 
terface [29j; the range of heights considered here allows 
us to neglect gravitational [lfl] and disjoining [13] pres- 
sures. Assuming that the height gradients are small, the 
lubrication approximation applies, and the local Laplace 
pressure is given by p(x,t) » —jd^h. Though the as- 
prepared stepped films have steep gradients, they also 
contain large curvature gradients in the same regions. 
These high-gradient regions then flow quickly compared 
to the experimental annealing times and the lubrication 
approximation is valid for all experiments presented here. 
As a stepped film levels, the bump represents a region of 
high pressure relative to the flat regions on the thick side 
of the film, while the dip has a low pressure relative to 
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flat regions on the thin side. Thus, there are pressure 
gradients along the region of transition between h\ and 
hi + h 2 and the flow continues. 

To understand the data of Fig. [Ha) quantitatively, we 
combine the lubrication approximation and the Stokes 
equation of viscous flows that leads to d x p = "qd^v, where 
v{x,z,t) is the horizontal velocity and z is the vertical 
coordinate 31-33J. Since the longest relaxation time of 
any polymer used here is of order 100 s, and all leveling 
times are significantly longer than this time, viscoelastic 
effects can safely be neglected. Assuming no slip at the 
solid-liquid interface and no stress at the liquid-air inter- 
face gives a Poiseuille velocity profile along z. Invoking 
conservation of volume and using the Laplace pressure 
introduced above leads to the thin film equation 



d t h 



^d x (h 3 d*h) = o 



(i) 



This differential equation can be nondimensionalized 
by letting H = h/h 2 , X = x/xq and T — jh^t/SijXQ, 
where Xg is a typical horizontal length scale of the prob- 
lem [31(. Furthermore, introducing the variable 



U = 



X 

fl/4 



1/4 



(2) 



one can show that self-similar solutions of the first kind 
of the form H(X,T) = F{U) exist [IS Using a 

numerical scheme (33, HI] , we have shown that this self- 
similarity is satisfied for times T > 10~ 4 for a profile 
with hi — h 2 . We have furthermore verified that the 
self-similarity is not sensitive to the initial height profile, 
provided that the heights at U = ±oo are constant and 
different. Therefore, if the interface profile is self-similar, 
a rescaling of the horizontal axis must collapse all of the 
experimental data. In Fig. &b) we show the same data 
as in Fig. 0(a) but plotted as a function of x/(h 2 t) 1 ' A . 
Having done so, we see that the data for all three molecu- 
lar weights do collapse onto three individual curves, thus 
demonstrating the self-similarity of the evolution through 
time. Examining the scaled curve for the 192 kg/mol PS 
in Fig. closely, one can see that the early time data 
do not collapse perfectly onto one single curve. This re- 
sult is due to the fact that it takes longer times to reach 
the self-similar solutions to Eq. (fTJ) for higher molecular 
weights, at a given temperature. 

In the previous rescaling, the horizontal length scale 
characterizing the transition between undisturbed film 
heights is different for each of the three data sets. Since 
the samples all started with the same initial condition, 
they have different leveling speeds because the viscosity 
varies with M w . 

To obtain precise measurements of these leveling 
speeds, we numerically solve the dimensionless form of 
Eq. (fTJ) with a stepped initial condition of aspect ratio 



r = h±/h 2 . We then fit the computed self-similar pro- 
file, F r (U), to the corresponding rescaled experimental 
profile. Figures Ela) and EJb) show the rescaled mea- 
sured profiles for two 15 kg/mol PS stepped films with 
different {hi, h 2 } values, as well as the corresponding fits 
to the numerical solutions of Eq. (fTJ) [Hj]. Despite the 
initial transient flow not described by the lubrication ap- 
proximation, all experimental profiles presented here ap- 
proach the self-similar solutions obtained from Eq. (flj. 
In the fitting procedure, the only free parameter is a hor- 
izontal stretch of the computed profile given by the fac- 
tor (7/377) x / 4 according to Eq. {2]). This procedure thus 
yields a measured value of 77/7. Using 7 « 30 mJ/m 2 [3§| 
for all polymers here, we get 77 = 7.3 x 10 3 , 1.4 x 10 5 , 
1.1 x 10 6 , and 9.1 x 10 6 Pas for M w = 15.5, 55.5, 118, 
and 192 kg/mol PS at 140 °C. These viscosities follow 
the expected molecular weight dependence of the poly- 
mer melt viscosity [40l | and are in agreement with viscosi- 
ties measured in the bulk [27]] ■ We have measured 77/7 
for nine additional 15.5 kg/mol PS stepped films all with 
different {hi,h 2 } values. We find that the value obtained 
is independent of {hi, h 2 }, and the ratio of the standard 
deviation to the mean is 0.12. 

Having exhibited the self-similarity of the profiles and 
having set up a robust method to measure the capillary 
velocity, we now turn to the study of the dissipation law 
in stepped films. The excess free energy per unit length 
along the dimension of invariance, AJ 7 , is dominated by 
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FIG. 3: (a) Rescaled measured profile for a 15 kg/mol PS 
stepped film with {/ii,/i2} = {30,193} nm after 10 min of 
annealing at 140 °C; (b) an identically annealed stepped film 
with {/ii,/i2} = {174,32} nm. The lines in (a) and (b) are 
horizontally stretched numerical solutions from Eq. {TJ using 
Refs. [37|, [1^1 . (c) Relative excess contour length of the data 
in Fig. [5] as a function of time; see Eq. 
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FIG. 4: Relative excess contour length as a function of time, 
the geometrical and physical parameters of the stepped films, 
and having made use of /3 = —0.42. Data points represent 
15 kg/mol films at t = 10 min with different {/ii,/i2} values 
(0.16 < r < 5.6), as well as the temporal series of Fig. [3jc). 
The solid line is the theoretical prediction of Eq. (JSJ with no 
free parameter (£ = 0.16). 



the interfacial contribution. For small slopes 



AT = 7 AL 



dx (d x hy 



(3) 



where AL is the excess contour length of the profile with 
respect to the flat limit at infinite time. In Fig. 02[c), we 
plot the relative excess contour length as a function of 
time for the three samples of Fig. [2j It is evident that 
the individual data sets obey a power law in time. Refer- 
ring also to the two stepped films shown in Figs.[3ja) and 
EJb), we see that the contour lengths of these two pro- 
files are not identical despite having been annealed under 
the same conditions. Therefore, there are temporal and 
geometrical dependencies on the energy dissipation. 

In the following, we develop a theoretical model in 
order to understand the results of Fig. [3l Making the 
quantities in the integrand of Eq. ([3]) dimensionless as 
presented above and invoking the self-similarity of the 
profile, we find that the relative excess contour length 
remaining after time t satisfies 



AZ = C(r)(^ 



with C{r) 



3V4 
2 



dU [Fl(U)] 



(4a) 
(4b) 



where C(r) is determined by the sample geometry 
through the dimensionless self-similar profile F r {U). 

The stepped film problem contains three relevant com- 
binations of heights: hi for the step height which pro- 
vides the typical driving force for leveling, as well as hi 



and h\ + hi for the thicknesses of the thin and thick re- 
gions of the film. Thus, having already extracted the hi 
dependence in Eq. (|4aj) . the simplest ansatz for a height 
dependence of C (r) is a power law of the form 



C(r) = C(l) r° 



1 



P 



(5) 



This scaling expression connects the case of arbitrary as- 
pect ratio, r, to the r = 1 geometry; normalization to 
r = 1 is ensured by the factor of 2. In order to determine 
the two exponents, we consider the limit r » 1, where 
the step is a small perturbation of the flat film. In this 
limit, Eq. flU can be linearized and solved analytically 
and we find C(r) w £ r~ 3 / 4 , where £ « 0.16 [4lJ. By 
comparison to Eq. §5§ when r 3> 1, we get a + f3 = —3/4 
and (3 = log[C(l)/£]/log(2). Then, using the numer- 
ical solutions of Eq. (JlJ presented above, we calculate 
C(l) w 0.12 and thus (3 -0.42. Finally, we obtain a 
master expression for the relative excess contour length 



AL 
~h7 



3/4 
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(6) 



Figure 0] shows the relative excess contour length for 
all samples and times presented in this study as a func- 
tion of the combination of {/ii, hi}, 7/77 and t suggested 
by Eq. ([5]). The solid line corresponds to the theoretical 
prediction of Eq. (|6]) with no free parameter. The agree- 
ment between the data and Eq. © validates the dissipa- 
tion law, as well as the general scaling to the r = 1 sym- 
metric case, within the range of our experimental data. 
As discussed above, the dissipation law is modified for 
t < 100 min for the 192 kg/mol PS, where self-similarity 
may not be valid. Although we have not accessed these 
limits, we expect deviations when r — > for which the 
dependency on hi may differ, and when t — ¥ for which 
the lubrication approximation is not valid. 

In conclusion, we have demonstrated that the thin film 
equation captures the time evolution of stepped polymer 
films for a wide range of aspect ratios and material prop- 
erties. With the use of a numerical solution of this equa- 
tion, it is now straightforward to precisely measure the 
capillary velocity of any nonvolatile fluid prepared as de- 
scribed. We further demonstrated the self-similarity of 
the evolution of the profiles. Finally, details of how the 
excess surface energy is dissipated by viscosity over time 
have been obtained. We have determined that the surface 
energy decreases with a — 1/4 power law in time. The rate 
at which energy is dissipated depends on the capillary ve- 
locity and on the heights of the initial stepped film. In 
particular, we have shown through a master expression 
that the evolution of any stepped film can be rescaled to 
the hi — hi symmetric step. 
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